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Global variables such as the transverse energy and the charged particle multiplicity and their ratio 
are evaluated, in a statistical model with expansion, for heavy-ion collisions from AGS to RHIC at 
' -y/sjvjv ~ 200 GeV. Full description of decays of hadron resonances is applied. The predictions of 

' the model done at the freeze-out parameters established independently from observed particle yields 

, and pt spectra agree well with the experimental data. However, some (explicable) overestimation 

■ of the ratio has been found. 

O ' PACS numbers: 25.75.-q, 25.75.Dw, 24.10.Pa, 24.10.Jv 
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I. INTRODUCTION 



> 

lO i The statistical model has been applied successfully in description of the soft part of particle production in heavy-ion 
O ■ collisions 0, 0, Namely, the particle yield ratios and pr spectra were fitted 

, accurately with the use of only four parameters (for particle ratio fits two parameters are enough: the temperature T 
^ ■ and the baryon number chemical potential hb)- Now, those parameters will be used to evaluate global observables: 
I the transverse energy density dEx/dri, the charged particle multiplicity density dNch/drj and their ratio (for details 
. see 53)- The advantage of such an approach is based on the fact that the transverse energy and charged particle mul- 
[ tiplicity measurements are independent of hadron spectroscopy (in particular, no particle identification is necessary), 
therefore they could be used as an additional test of the self-consistency of a statistical model. 
The experimentally measured transverse energy is defined as 



ET = Y.E,-sine,, (1) 



, where 6i is the polar angle, Ei denotes Ei — tun [rriN means the nucleon mass) for baryons and the total energy Ei 



for all other particles, and the sum is taken over all L emitted particles Additionally, in the case of RHIC at 
^/sNN = 200 GeV, Ei + rriN is taken instead of Ei for antibaryons IItII. 

The statistical model with single freeze-out is used (for details see [131 )• The model very well reproduces ratios and 
Pt spectra of particles measured at RHIC H, |^ . The main assumption of the model is the simultaneous occurrence 
of chemical and thermal freeze-outs, which means that the possible elastic interactions after the chemical freeze-out 
are neglected. The conditions for the freeze-out are expressed by values of two independent thermal parameters: T 
and ^B- The strangeness chemical potential fis is determined from the requirement that the overall strangeness of 
the gas equals zero. 

The actually detected (stable) particles have two sources: (a) a thermal gas and (b) secondaries produced by decays 
and sequential decays of primordial resonances. The distributions of particles from source (a) are given by a Bose- 
Einstein or a Fermi-Dirac distribution at the freeze-out point. The distributions of secondaries (source (&)) can be 
obtained from the elementary kinematics of a many-body decay or from the superposition of two or more such decays 
(for details see 0,^3)- ^^'^ following, all possible (2-, 3- and 4-body) decays with branching ratios not less than 
1% are considered. Also almost all cascades (with the exclusion of a very few like 4- or 5-step ones such that they 
proceed via at least one 3- or 4-body decay) are taken into account. It should be stressed that all contributions from 
weak decays are included in the presented analysis. 
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II. THE FOUNDATIONS OF THE SINGLE-FREEZE-OUT MODEL 

The foundations of the model are as follows: (a) the chemical and thermal freeze-outs take place simultaneously, 
(&) all confirmed resonances up to a mass of 2 GeV from the Particle Data Tables ^18^ are taken into account, (c) a 
freeze-out hypersurface is defined by the equation 



T = ^t'^ — r'^ — r"^ — r1 = const , (2) 
{d) the four- velocity of an element of the freeze-out hypersurface is proportional to its coordinate 



M - _ 



w 

(e) the transverse size is restricted by the condition r = y'r^ + < pmax- In this way one has two additional 

parameters of the model, r and Pmax, connected with the geometry of the freeze-out hypersurface. 
The maximum transverse-flow parameter (or the surface velocity) is given by 

\/l + {Pmax/rY 

The invariant distribution of the measured particles of species i has the form IToj 

dN, 



d^PT dy 



p'^da^ f,{p-u) , (5) 



where da^ is the normal vector on a freeze-out hypersurface, p ■ u — p^u^ , is the four-velocity of a fluid element 
and fi is the final momentum distribution of the particle in question. The final distribution means here that i s the 
sum of primordial and simple and sequential decay contributions to the particle distribution (for details see [131 )• 
The pseudorapidity density of particle species i is given by 

dNi [ 2 dy dNi _ [ ^ p dN, 



d-q J drj d'^pr dy J d^pr dy 

Analogously, the transverse energy pseudorapidity density for the same species can be written as 

d-q J p d-q d^pr dy J Ei d^pr dy ' 

For the quantities at midrapidity one has 



d-q 



, ^ Pt -75 1 , (o) 

mid J d'^pT ay -rriT 



dE' 



TA 



d-q 



J dpTPT " ' * = nudeon 



J d^pT PT d4%v ' ^ ^ nudeon 



(9) 



where Wc.m.s is the velocity of the center of mass of two colliding nuclei with respect to the laboratory frame (only for 
RHIC Vc.m.s = 0). Note that for RHIC at ^snn = 200 GeV there is the third possibility in Eq. 0: if i = a-ntinudeon, 
the analogous formula as for i = -nudeon but with mT + -m^ in the numerator holds. 
The overall charged particle and transverse energy densities can be expressed as 
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ch 
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(10) 



dET 
drj 



mid 
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drj 



mid 



(11) 



where A and B {B C A) denote sets of species of finally detected particles, namely the set of charged particles 
B = {tt+j tt^, , K^, p, p}, whereas A also includes photons, Kl, n and n 



III. RESULTS 

The general scheme reviewed in the previous section was formulated orig inally for RHIC 's','^,'!^ and then applied 
for SPS 01 . Here, this method will be used also for the AGS case. But the different model of the freeze-out 
hypersurface was applied for the description of pr spectra there 0,0- In that model (for details see the freeze- 

out happens instantaneously in the r direction, i.e. at a constant value oft {not at a constant value of r as here). The 
parameters connected with the expansion are the surface velocity /J™"^ and Pmax- Therefore, the implementation of 
values of obtained within that model into the presented one is entirely ad hoc, nevertheless it works surprisingly 
well. To put values of from Q, || into formulae of sect.m one should invert Eq. to obtain 



^ fOmax 

Hmax hJ± (121 

T ~ y/l - (/3™'^^)2 ■ 

It should be recalled here that the value of r itself is not necessary to calculate the transverse energy per charged 
particle, since this parameter cancels in the ratio. 

The final results of numerical estimates of dET/drj\mid and dNch/dr]\mid together with the corresponding exper- 
imental data are listed in Table|l| To make predictions for the AGS case it has been assumed that the maximal 
transverse size Pmax equals the average of radii of two colliding nuclei and the nucleus radius has been expressed as 
Ra = roAs , tq = 1.12 fm. Generally, the estimates agree well with the data. However, for RHIC the 11% — 16% 
underestimation of the charged particle density has been found. But this simply reflects the existing inconsistency 
in measurements of the charged particle multiplicity at RHIC. Namely, the sum of integrated charged hadron yields 
|20l |. after converting to dNch/drj 0|, is substantially less then the directly measured dNch/dr]\mid Hj. This is shown 
explicitly in the last column of TableQ] But values of the sum agree very well with the model predictions. Since the 
geometric parameters were established from the fits to the same px spectra, the agreement had to be obtained. Also 
for AGS the results agree qualitatively well with the data, in spite of the roughness of the method applied for this 
case. The overall error of evaluations of transverse energy and charged particle densities is about 0.5% and is caused 
by: (a) omission of the most complex cascades; (6) simplifications in numerical procedures for more involved cascades. 
The velocity of the center of mass of two colliding nuclei, Vc.m.s, equals: for RHIC, 0.994 for SPS Pb-Pb collisions 
at 158 • A GeV, 0.918 for AGS Au-Au collisions at 11 • A GeV and 0.678 for AGS Si-Pb coUisions at 14.6 • A GeV. 

Values of the ratio dEx / d'q\mid/ dNch/ dri\mid can be also calculated. They are collected in TablcHll together with 
the corresponding data. The overall overestimation of the order of 15% has been obtained. In the RHIC case this is the 
result of the underestimation of dNch/dr]\midi which has been explained earlier. But when in the denominator of the 
experimental ratio, dNch / drjlmid from the summing up of integrated hadron yields is put, the theoretical predictions 
agree very well with the data. Note that the similar inconsistency in charged particle measurements could have also 
been the origin of the discrepancy between model and experimental values of dNch/drj\mid seen in the AGS Si-Pb 
case. For SPS, the result agrees with the experimental value within errors. The overall error of model evaluations of 
the ratio is less than 1%. These results have been also depicted together with the data in Fig.Q] One can see that 
the relative positions of theoretical points agree very well with the data, they are shifted up only and this is the effect 
of the overestimation discussed earlier. 



IV. CONCLUSIONS 



The expanding thermal hadron gas model has been used to reproduce transverse energy and charged particle 
multiplicity pseudorapidity densities and their ratio measured at AGS, SPS and RHIC. The importance of the present 
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TABLE I: Values of dEr / dr]\mid and dNch / drj\mid calculated in the framework of the statistical model with expansion. In 
the first column thermal and geometric parameters are listed for the corresponding collisions. In the third and last column 
experimental data for the most central collisions are given. 



Collision case 


dET/dn\m^d [GeV] 


dNch/drtlmid 




Theory Experiment 


Theory Experiment 


Au-Au at RHIC at y'sjviv = 200 GeV: 
T = 165.6 MeV, hb = 28.5 MeV 
Pmax = 7.15 fm, T = 7.86 fm (/JJ'"" = 0.67) [11] 

Au-Au at RHIC at ^^snn = 130 GeV: 
T = 165 MeV, fiB = 41 MeV 
p^„, = 6.9 fm, T = 8.2 fm {PT"" = 0.64) [12] 


585 " 597 ± 34 [17] 
507 503 ± 25 [16] 


589 699 ± 46 [17] 
579 ± 29 ' 
[22] 

555 622 ± 41 [21] 

568 ± 47 

476 464li!^ [2^ 


Pb-Pb at SPS: 
T = 164 MeV, /ib = 234 MeV 
Pn.a. = 6.45 fm, r = 5.74 fm {PT"'' = 0.75) [13, JJ 


447 363 ± 91 [2^] 


Au-Au at AGS: 
T = 130 MeV, fiB = 540 MeV 

^rna. ^ q q^^^^ ^^^^ ^ g 53 fm [1 , J 


224 ^ 200 [24] 


271 ^ 270 [25] 


Si-Pb at AGS: 
T = 120 MeV, iiB = 540 MeV 
^max ^ Q_54^ ^^^^ ^ 5 Q2 fm [1, J 


57 ^ 62 [25] 


91 ^ 115 - 120 
[25] 


"For the modified definition of Ej-, «.e. Ei + mjv is taken instead of Ei for antibaryons, see Eq. tTl. 
'For the charged particle multipUcity expressed as the sum of integrated charged hadron yields. 




TABLE II: Values of the ratio dET / dr]\mid/dNch/drj\mid calculated in the framework of the statistical model with expansion. 
In the last column experimental data for the most central collisions are given. 


Collision case 


dET/dri\mid/dN^h 


/d-q\mid [GcV] 




Theory 


Experiment 


Au-Au at RHIC at ^snn = 200 GeV 
Au-Au at RHIC at ^Jsnn = 130 GeV 
Pb-Pb at SPS 


0.99 " 
0.91 
0.94 


0.87 ±0.06 [12 
1.03 ±0.08 

0.81 ±0.06 [16] 
0.89 ±0.09 ' 

0.78 ±0.21 [23| 


Au-Au at AGS 


0.83 


0.72 ±0.08 [25] 


Si-Pb at AGS 


0.63 


0.52-0.54 [2^ 



"For the modified definition of Ej-, «.e. Ei + mjv is taken instead of Ei for antibaryons, see Eq. tTl . 

''Author calculations with the use of experimental values given in Tableland the denominator expressed as the sum of integrated charged 
hadron yields. 
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FIG. 1: Values of the transverse energy per charged particle at midrapidity for the most central collisions. Black dots denote 
evaluations of the ratio in the framework of the present model (the second column of TablelTB . Also data points for AGS |25| 
(a circle for Au-Au and a vertical bar for Si-Pb), SPS '23] (triangle), RHIC at ^/sNN ~ 130 GeV 16] (square) and RHIC at 
= 200 GeV [H (star) are depicted. For RHIC, points with the sum of integrated charged hadron yields substituted for 
the denominator are also depicted (crosses). 



analysis originates from the fact that the transverse energy and the charged particle multiplicity are independent 
observables, so they can be used as new tools to verify the consistency of predictions of a statistical model for all 
colliders simultaneously. The predictions have been made at the previous estimates of thermal and geometric freeze- 
out parameters obtained from analyses of measured particle ratios and pT spectra at AGS ,2., JJ, SPS 0, 0| and 
RHIC 1^ 10. The overall good agreement, not only of the ratio but also absolute values of dEx/d-q \mid and 
dNch/dr]\mid, with the data has been achieved. And the observed discrepancies can be explained reasonably. This 
strongly supports the idea that the thermal expanding source is responsible for the soft part of the particle production 
in heavy-ion collisions. In addition, the description of various observables is consistent within one statistical model. 

There are more arguments in favour of the above statement. One could think that this analysis is a kind of an 
internal consistency check of various measurements. And such a check could be done even in an model- independent 
way simply by integrating spectra of stable particles. But it can not be done without any external input. First, 
transverse momentum spectra are measured in limited ranges, so very important low-px regions are not covered by 
the data. Therefore, to obtain integrated yields some extrapolations below and above the measured ranges are used. 
In fact these extrapolations are only analytical fits without any physical reasoning, but contributions from regions 
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covered by them account for 25% — 40% of the yield 20]. On the other hand, a calorimeter acts very effectively in 
the low-pT range, namely pions with pr < 0.35 GeV/c, kaons with pT < 0.64 GeV/c and protons and antiprotons 
with pt < 0.94 GeV/c are all captured [l^. Since the very accurate predictions for the transverse energy density 
at midrapidity have been obtained (see Tabled , the present analysis can be understood as an undirect proof that in 
these unmeasurable px regions spectra are also explicable by means of the thermal source with flow and decays. 

Moreover, it is impossible to check the consistency of the transverse energy data because not all stable hadron 
spectra are measured. This mainly concerns neutrons and K^. Also it is impossible to extract hadron decay photons 
from the photon data. And again, the very good agreement of model estimates of the transverse energy density at 
midrapidity with the data can be interpreted as the strong argument that the production of neutral stable particles 
can be described in terms of the expanding thermal source with superimposed decays. 

And as the last remark: in opposite to the transverse energy, there is some inconsistency (of the order of 10%) of 
the independent measurements of the charged particle multiplicity with the corresponding sums of integrated charged 
particle yields at RHIC (see sect. But at the present stage of investigation it is difficult to judge whether this 
inconsistency has the physical or experimental reason. 
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